UNCLASSIFIED 


_ AD  NUMBER _ 

AD913120 

LIMITATION  CHANGES 
TO: 

Approved  for  public  release;  distribution  is 
unlimited. 


FROM: 

Distribution  authorized  to  U.S.  Gov't,  agencies 
only;  Test  and  Evaluation;  17  SEP  1973.  Other 
requests  shall  be  referred  to  Electronic 
Systems  Division,  ATTN:  ESD-TR-2,  Hanscom  AFB, 
MA  01730. 


_ AUTHORITY 

ESD  ltr  dtd  4  Dec  1973 


THIS  PAGE  IS  UNCLASSIFIED 


Unclassified 


Project  Report 


C.  R.  Berndtson 
R.  H.  French 

D.  E.  Nessman 


Prepared  for  the  Advanced  rch  Projects  Agency, 

the  Department  of  the  Army,  and  the  Department  of  the  Air  force 
ir.der  Electronic  Systems  Division  Contract  F19628-70-C-0230  by 

Lincoln  Laboratory 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 

. 

Lexineton,  Massachusetts 


t 


Unclassified 


84 


MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
LINCOLN  LABORATORY 


DATA  REDUCTION  PROGRAM  DOCUMENTATION  ALTAP 


(EFFECTIVE:  APRIL  1971). 


\  \  /  \ 

/  /O)  C  RyBERNDTSON j  J 

Gi— r  oa 


r.  h.  French 
•  d.  e  Aiessman 

l  "it _ s 

Hi  tco-Ford  Corporation 
Editors 


/)1\  f/j  d  tr 


PROJECT  REPt«TjiPA-229-4((RSP) 


on 


2  AP RICMM71 


Distribution  limited  to  U.S.  Gov't,  agencies  onlyf 
Test  and  Evaluation;*  ^  SEfiJ973*  Other  iequests 
for  this~Hocument  must  be  referred  to 


with  specific  prior  approval  of  ESD-TR-2. 

L'  Q- 


LEXINGTON 


434 


W  07? 

Unclassified 

xoj  t 


MASSACHUSETTS 


Unclassified 


The  work  reported  >n  this  document  was  performed  at  Lincoln  Laboratory, 
a  center  for  research  operated  by  Massachusetts  Institute  of  Technology. 
This  work  was  sponsored  by  the  Advanced  Research  Projects  Agency  of  the 
Department  of  Defense  (ARPA  Order  600\  the  Department  of  the  Army,  and 
the  Department  of  the  Air  Force  under  Air  Force  Contract  F19628-70-C-0230. 


ii 

Unclassified 


FOREWORD 


This  is  the  fourth  rep^~t  in  the  Data  Reduction  Program  Documentation 
series.  It  is  dated  according  » the  date  of  completion  of  the  documentation.  No 
implication  is  made  that  this  program  will  not  subsequently  be  modified,  amended,  or 
superseded;  on  the  contrary,  the  history  of  radar  data  processing  is  one  of  continuous 
evolution  of  techniques,  and  it  is  unrealistic  to  assume  that  steady -state  has  been 
reached.  The  PA -229  series  is  being  published  for  the  convenience  of  interested 
parties,  and  Lincoln  assumes  no  responsibility  for  the  correctness  of  the  information 
presented,  nor  for  its  currency. 

The  preparation  of  reports  in  this  series  is  under  the  Editorship  of  Charles  R. 
Bemdtson  of  Lincoln,  and  of  D.  Nessman  and  R.  French  of  Philco-Ford  Corporation. 
Inquiries,  suggestions,  corrections,  criticisms,  and  requests  for  additional  copies 
should  be  directed  to  C.  R.  Bemdtson. 

The  principal  contributor  to  this  report  was  A,  J.  Poirier  (Philco-Ford). 

Due  to  the  intricate,  evolutionary  manner  in  which  the  programs  came  into  being, 
the  editors  regret  that  it  is  in  general  impossible  to  give  due  credit  to  all  -- 
mathematicians  or  radar  analysts  or  programmers  --  who  contributed  to  the 
definition  and  writing  of  the  programs. 


Alan  A.  Grofrietstein 
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COMMON  SYMBOLS  AND  ABBREVIATIONS 

(The  units  given  for  certain  quantities  are  the  units  commonly  used  for  those 
quantities,  unless  otherwise  noted. ) 


ADT 

Alt 

APS 

ARS 

Avg 

Az 

CADJ 
C -band 


ALCOR  Data  Tape 
Altitude  (km) 

Average  Pulse  Shape 
ALT  AIR  Recording  System 
Average,  Averaging 
Azimuth  (deg) 

Adjusted  Calibration  Constant  (db) 

ALCOR  frequency,  5664  MHz  (NB)  and  5667  MHz  (WB) 

Elevation  (deg) 

End  of  File 


GMT 


Greenwich  Mean  Time 


h  Hours 

Hz  Hertz 

in  Inches 

LC  Left  Circular  Polarization 

min  Minutes 

NB  Narrow  Band 

NRTPOD  Non-real  Time  Precision  Orbit  Determination  Program 


POD 

Phase 

PRF 

PRI 

pps 

pts 


Project  PRESS  Operation  and  Data  Summary  Report 

Presented  in  deg 

Pulse  Repetition  Frequency  (pps) 

Pulse  Repetition  Interval  (s) 

Pulses  per  second 
Points 

Range  (km) 

Range  Rate  (km/s) 

Radians 

Right  Circular  Polarization 
Radar  Cross  Section  (dbsm) 


vi 


s  Seconds 

SD  Standard  Deviation  of  Wake  Velocity 

w 

T  Time 

TAL  Time  After  Launch  (s) 

UHF  ALT  AIR  Frequency;  415  MHz 

V  Velocity 

V  Doppler  Velocity 
d 

V  Mean  Wake  Velocity 
w 

VHF  ALTAIR  Frequency;  155.5  MHz 

WB  Wide  Band 

0  Total  Off-a^iS  Angle  (deg) 

A.  Wavelength 

*  Denotes  Multiplication 
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F i.v W  DIAGRAM  SYMBOLS 


PROCESS,  ANNOTATION 


DECISION 


TERMINATOR 


SUBROUTINE:  where  NAME  is  the  entry 
cell  Into  the  subroutine 


CONNECTOR:  where  P  specifies  a  page  In  the 

flow  diagram,  and  L  designates 
a  statement  number  In  the  'ogram 
listing  or  a  reference  point  In  the 
flow  diagram 

CONNECTOR:  where  X  Implies  a  continuation 
of  the  diagram  to  the  next  page 


INPUT/OUTPUT  OPERATION 


MAGNETIC  TAPE 


PUNCHED  CARD 


DISK 
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ALTAP 


I. 


PURPOSE  AND  UTILIZATION 


A. 


B. 


C. 


Source  of  Data 


ALT  AIR 


Data  Input 


ALTAIR  transcription  tape 
pe  script  ion 


ALTAP  ca Igulate ^-ROSVtf  waking  targets  during  re-entry.  It  computes 
'or  up  to  120  range  gates.  It  is  normally  run  every  0. 1  s  with  the  data  averaged 


'over  0. 05  s. 


o.K 


\Output 

1.  Plots  of  RCS  vs  relative  range  at  altitude  increments  of  approximately 
2.5  km  or  for  every  averaging  interval.  Such  plots  are  non -coherent  average  pulse  shapes 
(APS)  which  show  the  position  of  the  wake  relative  to  the  body,  and  display  the  body  and 
wake  RCS. 

2.  Punched  cards  containing  TAL,  Alt,  gate  no.,  and  RCS  if  APS 
plots  were  obtained  every  2.5  km.*  These  cards  may  be  edited  and  modified  using 
results  of  coherent  data  analysis  programs  and  then  used  in  a  plotting  program  to 
produce  the  coherent  APS  shown  in  the  POD. 

3.  Punched  cards  containing  TAL,  Alt,  and  RCS  of  the  peak  gate. 

These  cards  may  be  edited  and  modified  using  results  of  coherent  data  processing 
programs  and  then  used  in  a  plotting  program  to  produce  a  peak  wake  plot. 


If  APS  plots  were  obtained  every  averaging  interval,  no  punched  cards  are  produced. 


'hjla&cJ l 


1 


120  gates. 


4. 


Listing  of  RCS  every  averaging  interval  for  a  maximum  of 


5.  Listing  of  RCS  i 


in  addition  to  dbsm  can  be  requested. 


* 


t 


II. 


DESCRIPTION 


ALTAP  computes  average  RCS  for  selected  gates  and  time  intervals.  ALTAP 
will  process  only  one  data  channel  on  one  target  per  run.  The  averaging  interval  may 
be  any  value  larger  than  the  PRI.  The  program  averages  in  m2  and  then  converts  to 
dbsm  for  printouts  and  plots. 

The  altitude  on  the  transcription  tape  has  a  resolution  of  1  km.  Therefore, 
the  altitude  used  in  ALTAP  is  computed  assuming  a  spherical  earth  as  follows: 

Alt  =  (R2  +  Re2  +  2RResn  El)*  -  Re 
where  Re  =  radius  of  earth  (6378. 145  km). 

A  number  of  input  parameters  and  transcription  tape  parameters  are  checked 
for  validity  before  processing. 

The  main  program  checks  the  following  input  parameters: 

0  <  I  POL  <  4 
IPAT  =  1  or  2 
TAVG  *  0 
INTARG  *  0 
NRG  *  0 

2 

Subroutine  ALREAD  makes  a  number  of  other  checks  on  transcription  tape 
parameters  which  are  described  in  detail  in  the  ALREAD  description.  For  some 
errors  (missing  format  tables;  end  of  file;  target  no.,  sampling  pattern,  or  polarization 
not  on  tape)  information  is  returned  to  main  program  for  decision  to  terminate. 
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III. 


OPERATION 


A.  Input 

Start  and  stop  times  (GMT) 

Averaging  interval  and  skip  time* 

Target  and  sampling  pattern  numbers 

Specified  set  of  range  gates 

Initial  gate  for  computing  peak  wake  RCS 

Alt  to  start  APS  plots 

Options  for  punching  and  plotting  data 

A  sample  ALTAP  input  is  shown  in  Appendix  A. 


CARD  1  (15A4) 

(Col.) 

1-60  TITLE  60  character  title  for  printout  and  plots. 

CARD  2  (213,  F7.3,  213,  F7.3,  4X,  515,  2F10.3) 

1-3  IH1  (13)  ) 

4-  6  IM1  (13)  >  Start  time  (GMT)  in  h,  min,  and  s 

7-13  ZSEC1  (F7.  3)\ 

14-16  IH2  (13)  ) 

17-19  IM2  (13)  >  Stop  time  (GMT)  in  h,  min,  and  s 

20-26  ZSEC2  (F7.  3)j 

31-35  NRG  Number  of  range  gates  (15) 

36-40  INTARG  Target  no.  (15) 

41-45  IPAT**  Sampling  pattern  in  which  initial  gate  is  located.  (15) 

t 

46-50  IPOL  Data  channel-  1  =  LC;  2  =  RC;  3  =  Az  error  ; 

4  =  El  error  t  (15) 
tt 

51-55  ISG  Location  within  IPAT  of  initial  gate  (15) 


*  Skip  time  is  the  time  in  seconds  from  the  end  of  one  averaging  interval  to  the  start 
of  the  next. 

**  Also  called  INPAT. 

1  VHF  transcriptions  only, 
tt  Also  called  ING  and  ISTGAT. 
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(Col.) 

5b -65 
66-75 

CARD 

1-  5 

6-15 

16-25 

26-35 

36-40 

41-45 

46-50 


TAVG  Averaging  interval  in  seconds  (F10.3) 

TSKIP  Skip  time  in  seconds  (F10.3) 

3(15,  3F10.3,  315) 


ISTGT  The  gate  location,  relative  to  ISG,  in  which  to  start  looking  for 

peak  wake.  If  0,  the  gate  location  is  set  to  1.  (15) 

REQAL  The  altitude  at  which  to  start  the  APS  plotted  at  altitude  intervals. 

(F10. 3)  This  must  be  on  tape  within  time  interval  being  processed. 

YMIN  The  minimum  ordinate  of  the  APS  plots  in  dbsm.  If  0,  the 

ordinate  is  set  to  -60.0  dbsm.  (F10.3) 

YMAX  The  maximum  ordinate  of  the  APS  plots  in  dbsm.  If  0,  the 

ordinate  is  set  to  +40.0  dbsm.  (F10.3) 

2 

ISQM  1  :  printed  output  in  m  and  dbsm. 

0  :  printed  output  in  dbsm  only.  (15) 

IOPT  If  >  0,  program  will  plot  every  averaging  interval.  (15) 

If  £  0,  program  will  plot  every  altitude  interval. 

IPUN  If  /  0,  program  will  punch  peak  wake  data.  (15) 

If  0,  punched  peak  wake  cards  not  produced. 


13.  Output 


LISTING 


GMT  and  TAL,  Alt,  and  R  for  mid -point  of  averaging  interval 
Frequency  and  polarization 

2 

Average  RCS  for  each  gate;  Average  RCS  (m  )  for  each  gate 
(when  requested). 

Number  of  pulses  used  in  averaging  interval  and  CADJ 

APS  PLOTS 

RCS  vs  relative  range  (m) 

PUNCHED  PEAK  WAKE  DATA 
Alt  (F10. 3),  RCS  (F10.3),  gate  no.  (15) 
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PUNCHED  APS  DATA  AT  ALTITUDE  INTERVALS  * 

Alt  (F10.3),  TAL  (F10.3) 

RCS  (F10.  3),  gate  no.  (15) 

(1  card  for  each  gate  selected) 

Sample  ALTAP  outputs  are  given  in  Appendix  B. 


*  Produced  only  when  10 PT  ^  0  and  REQAL  is  input. 


IV. 


PROGRAM  LIMITATIONS 


Start  time 

Must  be  on  tape 

Stop  time 

Must  be  on  tape 

NRG 

^  i 20  gates 

TAVG 

Must  be  larger  than  the  PRI 

TSKIP 

Cannot  be  negative 

INTARG 

Must  be  on  tape  within  scart  and  stop  times 

Length  of  run 

With  no  punched  cards:  no  limit. 

With  punched  cards:  £  300  averaging  intervals 

V. 


PROGRAMMING 


A.  TAPOS  (see  Appendices  C  and  D. ) 

TAPOS  is  the  control  section  of  ALTAP.  TAPOS  reads  the  input  cards, 
calls  ALREAD,  and  averages  the  data  returned.  TAPOS  also  calls  the  plot  routines, 
and  prints  and  punches  all  data. 

B.  PEAK  (see  Appendices  E  and  F.) 

PEAK  searches  an  array  of  RCS  to  find  the  largest  value;  it 
saves  this  value  and  the  gate  number  of  the  value.  PEAK  then  passes  these  values 
back  to  TAPOS.  A  location  within  the  array  of  RCS  may  be  specified  to  start 
the  search. 


ISPOT). 


The  call  statement  is  PEAK  (AVGAL,  NRG,  ISTGT,  IRGA,  AVGSX, 

INPUT 

AVGAL  Alt  * 

NRG  Number  of  range  gates 

ISTGT  Start  gate  for  peak  relative  to  ISG 

IRGA  Array  of  gate  nos. 

AVGSX  Array  of  RCS  for  gates 

OUTPUT 

ISPOT  Peak  gate  identification 

STORED  IN  COMMON 

INN  Running  number  of  averaging  intervals 

GTMAX  Array  of  peak  RCS* 

ALT  Array  of  altitudes* 

IGAT  Array  of  peak  gate  nos.* 


*  One  value  for  each  averaging  interval. 
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C.  ALTIT  (see  Appendices  G  and  H.) 

ALTIT  computes  Alt  by  using  R,  El,  and  the  radius  of  the  earth,  and 
returns  this  value  to  TAPOS. 

The  call  statement  is  ALTIT  (AVGAL,  AVGRG,  AVGEL). 

INPUT 

AVGRG  R* 

AVGEL  El* 

OUTPUT 

AVGAL  Alt* 

D.  REFC  (see  Appendix  J.  ) 

ALTIT  calls  REFC.  The  tropospheric  refraction  correction  sub- 

3 

routine,  REFC,  is  based  on  tropospheric  refraction  tables  in  PPP-36.  A  modified 
version  of  this  subroutine  is  now  in  use. 

The  call  statement  is  REFC  (E,  R,  DEE,  DRR). 

E  =  Uncorrected  El  (must  be  between  0  and  90  ) 

R  -  Uncorrected  R  (ft) 

DEE  =  El  tropospheric  correction 
DRR  =  R  tropospheric  correction  (ft) 

The  corrected  values  to  be  computed  after  exiting  >m  tl  REFC  subroutine  are: 

El  =  E-DEE 
R  (ft)  =  R-DRR 


*  For  midpoint  of  averaging  interval 


E. 


BLOTO 


BLOTO  plots  RCS  vs  relative  range  (m)  every  altitude  or  every 
averaging  interval.  The  ordinate  is  variable  through  the  optional  input  of  YMIN 
and  YMAX.  Nominal  values  for  these  are  -60.0  dbsm  and  +40.0  dbsm,  respectively. 

F.  TSPLIT  (see  Appendix  K.  ) 

TSPLIT  is  used  to  convert  time  from  GMT  total  seconds  to  h,  min,  s, 
and  decimal  fractions  of  s. 

The  call  statement  is  TSPLIT  (AVGTM,  IHM,  TRUN). 

INPUT 

AVGTM  GMT  total  seconds 

OUTPUT 

IHM  (1)  Hours 

IHM  (2)  Minutes 

TRUN  Seconds  and  decimal  fractions  of  seconds 

G.  ALREAD2 

ALREAD  is  the  Fortran  driver  for  the  machine  language  tape  reading 

routines. 

The  call  statement  is  ALREAD  (TSTART,  TSTOP,  TLIFT,  INTARG, 
INPAT,  I  POL,  NOPHA,  NPTS,  DFPG,  NEWPAS,  NRG,  ISTGAT). 

INPUT 

TSTART  Start  time  of  processing  (GMT  total  seconds) 

TSTOP  End  time  of  processing  (GMT  total  seconds) 

Target  number  to  be  processed 


INTARG 


Sampling  pattern  in  which  initial  gate  is  located 

Number  of  range  gates  to  be  processed 

Location  within  INPAT  of  initial  gate  wanted 

1  (only  RCS  data  wanted) 

Data  channel:  1  =  LC;  2  =  RC;  3  =  Az  error; 

4  =  El  error 

INPUT  and  OUTPUT  PARAMETERS 

NPTS*  Output:  number  of  pulses  of  data  returned 

Input:  must  be  initialized  by  calling  program 

before  each  call  to  ALREAD 

NEWPAS**  Cycle  and  error  pointer 

OUTPUT 

TLIFT  Lift-off  time  (GMT  total  seconds) 

DFPG  Frequency  and  polarization  (e.g.  VHF  LC) 

VALUES  STORED  IN  COMMON 

TIMES  Pulse  times  (GMT  total  seconds) 

XSPHA  RCS  and  phase  for  each  pulse  and  gate 

RANGKM  R 

ELSAV  El  (rad) 

CALOUT^  Adjusted  calibration  constant 

IRGA  Range  gate  array  associated  with  XSPHA 

NFPG  Frequency  code:  1  =  VHF;  2  =  UHF 


*Set  to  zero  for  first  call.  Set  to  number  of  saved  points  for  subsequent  calls. 

**Also  called  IAGAIN. 

^Also  called  CADJ. 


INPAT 

NRG 

ISTGAT 

NOPHA 

IPOL 
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H. 


REW 


4 

REW  is  an  entry  to  subroutine  BREADS  used  to  rewind  the  tape. 
J.  Plotting  System  Subroutines 

They  are  REREAD,  STOIDV,  and  PLTND. 
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APPENDIX  A 
ALTAP  INPUT 


'.I  Tc'iP  ri.1~.ii 


CARD  1 


0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  C  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0 

I  I  3  4  5  6  7  8  9  10  II  12  13  14  15  16  17  16  19  20  21  22  23  24  35  26  27  2(1  29  30  31  32  33  34  35  36  37  38  3S  40  41  42  43  44  4S  48  47  48  49  86  $1  52  S3  54  S5  56  37  56  59  GO  81  62  63  64  65  66  67  88  69  70  71  72  73  74  75  78  77  78  73  00 

i"'  nr  rn,r^  i?  nr  *>i .  no  ?r  19  1.  t  1  n.  ns  0.  nr, 

00  00  CARD  2 


000,00000000.00000,  .00000000000000000000000000000000000 


,uuu,uuuuuuu  l'  ,UUUUU,,UUUUUUUUUUUUUUUUUUUUUUUUUIIUUUUUUUUU,U,  0000000,  0,0000000 
I  2  3  4  5  6  7  9  9  10  11  12  13  14  15  16  17  18  19  20  21  ZZ  "3  24  25  26  27  Z0  29  30  31  32  33  34  3  6  36  3  7  31  39  40  41  42  41  44  45  4  6  4  7  4  0  49  50  51  52  53  54  55  SB  57  51  59  60  61  62  63  64  65  66  67  60  69  70  71  72  73  74  75  76  77  78  79  80 

^  "fin’n  cin' 0  '  0  1  CARD  3 


0000000000000,000000,0,0000000,0,00000000000,00000000000000000000000000000000000 
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APPENDIX  C 

TAPOS  PROGRAM  LISTING 


DIMENSION  DFPG (2) ,S0HSX  (120)  ,AVGSX ( 120)  ,  IDENT  ( 1 5)  ,IHM  (2)  , 
1GTMAX(300C)  ,ALT(3000)  ,IGAT(3000)  ,0PiOG{120)  ,ISTAR(121) 

EQUIVALENCE (I HM(1)  ,IHR)  ,  (I  HR  (2)  ,IBINJ 

COMMCN/RDCOMT/TIHES  (300)  ,  1SPHA  { 1  20 , 300)  ,  HANGAR  (300)  ,EIS£V(300)  , 
1CALOUT  (30  0)  ,IRGA  (120)  ,  NFFG 
COM  NON/PEEK/INN, GTMAX, AIT, IG AT 

DOUBLE  PRECISION  A VGTB ,S EC,S0MTH ,T  1  ,T2, TIRES, TOTI M ,TSI ART, TS70P , 
1TSV,  ZSEC,ZSEC1,ZSEC2,TI1FT 
DATA  IAST/'*'/ 

DATA  ISTA  R/1 2 1  * '  •/ 

DATA  IBLANK/'  '/ 

DATA  UPLOG/ 120*0.0/ 

DATA  SUMS  X/1 20*0 . 0/ 

DATA  AVGSX/120*0.0/ 

TCTIM (IH, IM, SEC) =E FLOAT (60* (60*1 H* IM) ) *SEC 
TST  OP  =  0 . 0 
IAGA IN=  0 
41  SUMTM=0.0 
I  SPOT =121 
NPTS=0 
INN=  0 

READ (5, 2C, END=901) IDENT 
20  FC  RM  AT  ( 1  5  A4) 

READ (5,60,END=901)IH1, IB1,ZSEC1 ,IR 2  ,1 B2f ZSEC2 , N RG, I'  T ARG f IPAT , IDO L 
1 ,ISQ,TAVG,TSKIP 

6  0  FORMAT  (2  (213,  F7.  3)  , 4X , 575 , 2F 10. 3) 

RFAD  (5,6  1)  ISTGT,RE0AL,THIN,YHAX,ISQ(1,  IOPT,IPUN 
61  FORMAT (I5,3F10.3,3I5) 

IF (ISTGT. EC.O) ISTGT* 1 
WFITE  (6,64) 

64  FORMAT (5X, 'THESE  ARE  YOUR  INPUT  CARDS') 

WPITE (6,62) IH1,IM1,ZSEC1,IH2,IM2,ZSEC2,NBG,INTARG, IPAT,IPOL, 

1 ISG,TAVG,TSKIP 

6  2  FORMAT (2  (2I3.F7.3)  ,4X,5I5,2M0,  3) 

WRITE (6,63) ISTGT,RECAL,YHIN,YMAX,ISQM,IOPT,IPUN 
63  FORMAT (I5,3F10.3,3I5) 

IF (NRG. EQ.O) GO  TO  550 
IF(INTARG.EQ.O) GO  TO  560 
IF ( (IPAT. EQ.O) .OR. (IPAT. GT. 3))  GO  TO  570 
IF ( (IPOL . EQ , 0) .OR. (IP0L.GT.4))G0  TO  5B0 
IF(TAVG.EQ.O.O) GO  TO  590 
70  TSTART=TOTIM (IH1,IH1,ZSEC1) 

IF ( (TSTART.GT.TSTOP) . AND. (IAOAIN. NE.44) ) GO  TO  72 
CALL  REW 
TAG AIN=  1 

72  TSTOP=TOTIM  (IH2,IM2,ZSEC2) 

IF ( (REQAL.EQ.0.0) . AND. (IOFT.EQ.O) ) GO  TO  75 
CALL  STOIDV (IDENT, 59,0) 

75  CALL  REREAD  (99,530) 

CCrJNT  =  0 
KOUNT=0 
INT-0 
T  1=T  STA  RT 
T2=T1*T AVG 

100  CALL  ALREAD(TSTART,ISTOP,TLIFT,INTARG,IPAT,IPOL,1,NPTS,DFPG, IAGA IN 
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1fNRG,ISG) 

If (I AGAIN. EQ. 44) GO  TO  510 
IF ( I AG A  IN. EQ. 55)  GO  TO  901 
IF(NPTS.EQ.O) GO  TO  91 
IF(IKT.EQ.O) HBITE  (6, 190)  IDENT,  INTARG 
190  FORMAT  ('  1  '  #30X'  ALTAIP  TAP  VERSION  19  FEB  1971  V/'MX 
115A9/31X, 'TARGET  NUMBER  «  ',15//) 

I  NT=  1 

DO  220  1=  1 f  NPTS 

1^0  IF(TIHES (I) .GT.T2) GO  TO  290 
IF  (T1.GT. TIMES (I)  )  GO  TO  220 
INPT=I 

IF(SUHTH.NE.O.O) GO  TO  180 
SUHTB=TIHES  (I) 

SUHRG=H ANGKH (I) 

SUHEL=ELSAV (I) 

CALSV=CALOOT (I) 

180  DO  2C0  K*  1 , NRG 

EXTEN=  (XSPHA  (K,I) /10.) 

IF(EXTEN.GT.75.0)GO  TO  295 
XSPHA  (KfI)  =  10. **EXTEN 
SUHSX  (K)  =  SUNS  X  (K)  ♦XSPHA(K,I) 

200  CONTINUE 

TSV=TIHES  (I) 

RSV=RANGKN (I) 

ESV=ELSaV  (I) 

CGUNT=CCUNT*1 
220  CONTINUE 

IF(IAGAIN.EQ.O)GO  TO  290 

NPTS=0 

GO  TO  100 

290  IF  (COUNT. NE. 0.0) GO  TO  280 
295  WRITE  (6 f 260) T2f EXTEN 

260  FORMAT  (/25X'AT  TIME  =  » F 1 2. 9 ,2 X • TH ERE  IS  A  TINE  GAP  CR  PAD  DA"1 
1  FXTEN  *  ' f  P10. 9) 

GO  TO  990 

280  DO  390  J=1,NRG 

UPLOG (J) =SUHSX  (J) /CCU NT 
IF (UFLOG (J) .GT.O. 0) GO  TO  300 
AVGSX (J) =99.99 
GO  TO  320 

300  AVGSX  (J)  =  10. *ALOG 10  (UPLOG  (J) ) 

320  SUHSX  (J)  =0.0 
390  CONTINUE 

AVGTH= (SUHTH+TSV) /2. 

AVGRG=(SUHRG*RSV)/2. 

AVGEl= (SUHEL+ESV) /2. 

IK=IFIX (COUNT) 

IKT=IABS (KOUNT-IK) 

IF( (MOUNT. NE.O) .AND. (IKT.GT.2) ) GO  TO  929 
KCUNT=IFIX (COUNT) 

KOUNT=COUNT 

CALL  TSPLIT (AVGTH#IHM, ZSEC) 

AVGTM=AVGTH-TLIFT 

CALL  ALTIT(AVGAL,AVGRG, AVGEL) 
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IF ( (REQAL.EQ.O.O) .AID. (IOPT. EQ. 0) ) GO  TO  358 

CALL  BLOTO(A7GAL,REQAL,RRGvYHI)f, THAI, COUNT, D PPG, IBGA, A VGSX , ZSEC, 
1IMIN,IHR, NFPG, A7GTM, IOPT) 

353  IF  (ISTGT.EQ.O)  GO  TO  359 

CALL  FEAK  (A7GAL,HBG,ISTGT, IRGA,A7GSX, ISPOT) 

359  WRITE  (6,3  60)  IHR,IMIN,ZSEC,A7GTM,  A7GRG,A7GAL,KOUNT,CALS V 

360  FORM  AT (/2X* TIME (GBT)  *  • ,213 , F8 . 4 ,2X , FI  2.4 , 2X, • R ANGE (KM)  =  •, 

1F10. 3,2X, 'ALT (KM)  =  * , F10. 3 ,2X, 14, •  PULSES ', 2X, • CADJ  =  *,F10.3) 

I  N=  1 

ISTAR  (ISPOT) =1 AST 
380  IOUT=IN* 1 7 

IF (IOUT.GT. NRG) IOUT=NRG 

WFITE  (6 , 4  CO)  (IRGA(L)  ,ISTAR(L)  ,L=IN,IGDT) 

400  FORM  IT ( IX 'RANGE  GATES • 3 X , 1 8  (15 , A 1) ) 

WRITE  (6 ,4  20)  DFPG,  (AVGSX(L)  ,L*IN,IOUT) 

420  FO RB  AT  ( 1 X , •  PPG  =  •  ,2A4 , 18F6.  1) 

IN= IOUT* 1 

IF  (IN. LE. RRG) GO  TO  380 
ISTAR  ( I S  P  OT )  =  I BL  A  H  K 
IF  (ISQH.LT. 1 ) GO  TO  440 

WRITE (6,421) IHR,IHIN,ZSEC,A7GTM, A7G RG , A7G AL , K  On NT 
421  FORMAT (/2X*  TIME  (GBT)  =  • , 21 3, P8 . 4 , 2X » TIM E (TS FC)  =  »,F12.4,?X, 

1 • RANGE  (KM)  =  • , F  10. 3, 2X , • A  IT  (KM)  =  •, F 10. 3, 2X, 14 , •  PULSES') 

NT=1 

427  ITUC=NI-»8 

IF (ITUO.GT.NRG) ITUO=NRG 
WRITE (6,422) (IRGA (L) , L*NI,ITUO) 

422  FORM  AT (IX' RANGE  GATES'9112) 

WRITE  (6,423)  (UPLOG  (L)  ,L*NI,ITUO) 

423  FORM  AT.  (IX* SQUARE  METERS  •  ,  1 P9E1 2 . 3) 

NI=ITUO+ 1 

IF  (N I. LE . NRG) GO  TO  427 
GO  TO  440 

424  KCU NT=0 

WRITE (6,425) AVGTM,IK 

425  FORMAT  (/2X‘ AT  TIME  =',F10.3,»  THERE  WAS  A  PR  F  CHANGE  THE  NUMBER  OF 
1  PULSES  WAS  • , 15) 

440  COUNT=0 

T 1=T2*TSKIP 
T2=T1  +  T  A VG 
SUMTM=0.0 
SUMRG=0. 0 
SUMAL=0.0 

IF(T2.LE.TIMES(NPTS))GO  TO  160 
IF (T2.GT.1STOP) GO  TO  510 
DO  460  K=INET, NPTS 
KNPT=K 

IF (II. LE. TIMES  (K) )  GO  TO  480 
460  CONTINUE 
480  ND  =  NPTS-KNPTH 

CO  5C0  N= 1 , ND 
H=KNPT*N-1 
TIMES (N ) =TI MES  (NL) 

RANGKH(N) =R ANGKM (NL) 

ELS  AV (N) =ELS A V (NL) 
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CAI.OOT  (N)  =C ALOUT  (NL) 

DO  5C0  L=1,NFG 
XSPHA  (L,N)=XSPHA(L,NL) 

SCO  CONTINUE 
NPT£=ND 

IF ( IAGAIN. NE. 0) 00  TO  100 
*>10  IF  (IPUN.EQ.O)  GO  TO  900 

WPITF  (9,511)  (ALT  (M)  ,GTHAX(M)  ,IGAT  (M)  ,(1*1,  INN) 

511  FOFMAT(2F10.3,I5) 

GO  TO  900 
550  WRITE  (6 ,555) 

555  FORMAT  (//2X' A  ZERO  VALUE. WAS  INPUT  FOR  NFG  THIS  IS  A  NO  NO') 

GO  TO  900 
560  WRITE(6,565) 

565  FORMAT (//2X' A  ZERO  VALUE  WAS  INPUT  FOP  THE  TARGET  »,  THE  DATA  EDIT 
1 OF  HAS  GOOFFED  AGAIN') 

GO  TO  900 

570  WRITE  (6 ,575) IPAT 

575  FORMAT  (//2X ' A  VALUE  OF',l5,'  WAS  INPUT  FOR  IPAT  THE  ONLY  LEGAL  VA 
HUES  FOE  IPAT  ARE  1,2,3') 

GO  T3  900 

550  WRITE  (6,585)  IPOL 

565  FORMAT (//2X'A  VALUE  OF',I5,»  WAS  INPUT  FOR  IPOL  THE  ONLY  LEGAL  VA 
1LUE3  FOR  IPCL  ARE  1,2, 3,9') 

GO  TO  900 
SOQ  WRTTE  (6,595) 

59c  FORMAT  (//2X ' A  ZERO  VALUE  FOR  TINC  CAN  NOT  WORK  IT  WILI  BE  SET  TO 
10.05  SECONDS  AND  THE  PRCGR  A  t  WILL  CONTINUE') 

TAVG=0. 05 
GO  TO  70 
90  C  IAG AIN=99 
GO  TO  41 

901  IF (  (REQAL. EQ. 0. 0) . AND.  (ICET. EQ. 0) )  GO  TO  902 
CALL  PLTND 
902  RETURN 
END 
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STOIDV 


initialtk 

SC4060  > 


COUK/T  =  0.Ot 

K0UNT  =  0 
T1  -  TSTART 
T2=  TS.TOP 


COUNT  IS  PULSES 
IN  CURRENT  INTERVAL. 
KOUNT  IS  PULSES 
TKi  PREVIOUS  INTERVAL. 
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APPENDIX  E 

SUBROUTINE  PEAK  PROGRAM  LISTING 


SUE  FOOT  I NE  PE AK  (  A  VGA! # NRG , ISTGT , IRG A , AVGSX , IS  POT) 
COMMCN/PEEK/INN,GT!1AX  ,  AIT,IGAT 

DIMENSION  IRGA  (NRG)  , AVGSX  { NRG)  ,GTM AX ( 3000) , ALT ( 3000)  , I GAT  (IP 00) 

INN=INN+ 1 

GMAX=-1000. 0 

DO  2 C  I -ISTGT, NRG 

IF  (  AVGSX  (I)  .  IE.  GMAX)  GO  TO  20 

GMAX^AVGSX  (I) 

IGAT  E  =  IPG  A  (I) 

ISP  01  =  I 
20  CONTINUE 

GTMAX (INN)  =  GMAX 
ALT (INN) = AVGAL 
IGAT  (INN)  =IGATE 
RETURN 
END 
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APPENDIX  F 

SUBROUTINE  PEAK  FLOW  DIAGRAM 


APPENDIX  G 

SUBROUTINE  ALTIT  PROGRAM  LISTING 


SUBROUTINE  ALT  IT  (  A  V GAL  ,  A  VGR  G  ,  A VG EL) 

DATA  BR,XKHFT,RE/. C  174 5 33 0003 048 , 6378 . 145/ 

RF= AVGRG/XRMFT 

AVGEE  =  A  VGF.L/DR 

CALL  REFC  (AVGEE,  RR  ,DEE,DRF) 

FANGE= ( RR-DRH) *XFMFT 
ELEV=AVGEL- ( (DEE-. 3) *  DR) 

ALT  =  SQRT  (RANGE* *2* PE* *2*  (2*  FaNGL’*PE*SI  N  (ELEV)  )  ) 

AVGAL=ALT-RE 

RETUFN 

END 
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APPENDIX  I 

SUBROUTINE  REFC  PROGRAM  LISTING 


S08  BODTINI  BZPC(F,H,DEI,DRR)  EFFECTIVE:  16  JUNE  WO 

DIHENSION  PE  (16,8)  ,  DB  <  16,8)  ,10(16)  ,10(8) 

DMA  DE/0.0  ,0.0  ,0.0  ,0.0  ,0.0  ,0.0  ,0.0  ,0.0 

10.0  ,0.0  ,0.0  ,0.0  ,0.0  ,0.0  ,0.0  ,0.0  ,0.0313,' 
20.0303,0.0292,0.0287,0.0282,0. 0272,0.0262,0.0253,0.0243, 0.0223, 
30.0214,0.0195,0.0171,0.0135,0. 0075,0.0  ,0.0937,0.0848,0.0770, 

40.0732,0. C694, 0. 0627, 0.C571 ,0.0522,0. 0480, 0. 04 1 2, 0. 03 85, 0. C337, 
50.0278,0.0205,0.0105,0.0  ,0. 185C,0. 1520,0.1250,0.1140,0.1050, 

60.0904,0.0795,0.0708,0.0636,0.0523  0.0478,0.0405, C. 03 23,0. 0229, 
70.01 14,0. 0  ,0.53 10,0.  3070,0.2120,0.  1830,0.1600,0.1280,0.  1060, 

80. 0 899, 0.C80, 0.06 12, 0.055 0,0. 04 55,0.  0354,0.02^,6,0.0120,0.0  , 

90.7550,0.  31 20,0.2400,0. 2020,0.  1750,0.  1370,0. 1 120 , 0. 0942,0.  081 1  , 
AO. 063 1,0.  0!  66,0.0466,0.0361,0.0350,0.0122,0.0  ,0.9120,0.4110, 

BO.  2560, 0.2  1 40,0.  18  40,0.  14 20 ,0.  1 1 50, P.  0967  0.083 0,0. 06 43,, 0.0575, 
CO. 0472,0. 0365,0. 025 2, 0. 0 122, D. 0  ,0.9700,0.4200,0.2600,9.2200, 

DO.  1900,0.  1460,0.  1 17  0,0.0980,0.084  0,0. 0653,0.0584,0.04  78,0.0  369, 
EO. 0254, 0.0123, 0.0  / 
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0, 

104. 

0, 

86. 

.6, 

73. 

9, 

57.1, 

46. 

4,  39. 

0, 

33.8, 

26. 

8. 

24. 

3, 

20. 

5 

A  16. 

6, 

13. 

0, 

9. 

8, 

8. 

4, 

446.0, 

172. 

0,105. 

0, 

87.4, 

74. 

0, 

58. 

0, 

46. 

6 

8  39. 

2, 

34. 

0, 

27. 

■  0, 

24. 

6, 

20.7, 

16. 

7,  13. 

0, 

10.0, 

8. 

4/ 

DATA 

ED, BTDEG/0. 

.01, 

2.0 

.0,5.0 

,6.0 

,8.0,10. 

0,12.0, 

,  14. 

0, 

18., 

20 

•  t 

124.  , 

30 

.  ,40 

•  1 

60.  , 

,90. 

.57 

.29578/ 

DATA 

RD/0. 

01 

,10. 

.,30 

I..60. 

,200., 

400. 

,1000. 

,2000./ 

IF  (R.  IE.  0.0)  GO  TO  300 
RG=R/6080 . 27 
DO  100  IED-2,15 
I* 1 7-IE  D 

IF  (E.GE.  ED  (I)  )  r,0  TO  120 
100  CONTINUE 
1=1 

120  DO  200  JRD=2,8 
J=10-JRD 

IF(RG,GE.  BD  (J)  )  GO  TO  220 
200  CONTINUE 
J=1 

220  IF( J.BQ. 8) GO  TO  340 

ZR=  ALOG  (BG/BD  (J)  )/ALOG(SD(J*1)/ID(J)) 

IP (E. IE. 0.0) GO  TO  320 

ZE=1L0G  (E/ED  (I)  )/ILCG  (ED  (1*1) /ED  (I)  ) 

DE1=  |(DE(I»1,J)-DE(I,J))*(1.-ZB)*(DB(I,J*1)-DE(I,J))*ZR)  *ZE 
DE2=  ( (DE  (I,J*1)-DE  (I,  J)  )  •  (1.-ZE)  ♦  (D!  (I*  1, 3*  1)  -DE  (I,J*  1) )  *ZE)  *ZR 
DFE=DE1 ♦DI2»DE (I, J) 

DB1  =  ((OB  (1*1, J)  -DB  (I,  J) )  *  (1.-ZI)  ♦  (DB  (I,  J*  1)  -DE  (I, J) )  *ZR)  *ZE 
DR2  =  ((Lft  (I,J*1)-DR  (I,J)  )•  (1.-ZI)  ♦  (DR  (1+  1,  J  +  1)  -DB  (I,  J*  1)  )  *ZE)  *ZR 
DRB=  (DR  14DB2+DR  (I ,  J)  ) 

GO  TO  400 


300  DIE=0.0 

DRH=0. 0 
GO  TO  400 

320  DIE  =  DE  (I,  J)  ♦  (DE  (I, J *1.  -DE  ( I,  J)  )  »ZB 
DBB=D?(  (I,  J)  »  (DB(I,  J*1)  -DR  (I,  J)  )  *BB 
GO  TO  400 

340  DELT* (E-ED(I))/(ED(I*1)-ED(I)) 

DEK=DEIT* (DE(I*1,J)-DE(I,J))  *DE(I,J) 
CRR  =  DELT*  (DR  (1*1,  J) -DR  (I,J)  )  *DB  (I,  J) 

400  RETOEN 
END 
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APPENDIX  K 

SUBROUTINE  TSPLIT  PROGRAM  LISTING 


SUBROUTINE  TSPLIT  (AVGTH,IHH,TRUN) 
DIHENSION  IHH  (2)  ,  DIVIDE  (2) 

DOUBLE  PRECISION  AVGTH,TBUN 
DATA  DIVIDE/3600. ,6C./ 

TRUN= AVGTH 
DO  20  1=1,2 
IHH  (I)  =TB UN/DIVIDE  (I) 
TFUN=TRUN-FLOAT (IHH  (I) ) *DIVIDE (I) 
20  CONTINUE 
RETURN 
END 


